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摘要 : 建立 了 以 无 数 个 无 限 深 势 阱 中 的 粒子 为 工 质 的 不 可 逆 量 子 斯 特 林 致 冷 循环 模型 。 势 阱 中 粒子 在 能 级 上 的 分 布 由 吉 布 斯 
分 布 函数 决定 。 该 致 冷 循环 由 两 个 等 能 量 过 程 与 两 个 等 势 阱 壁 宽 度 过 程 组 成 ， 具 有 高 低温 热源 间 的 热 漏 和 不 完全 回 热 因素 。 
导出 了 性 能 系数 、 致 冷 率 、 灶 产 率 和 生态 学 目标 函数 性 能 解析 式 ， 制 冷 率 与 生态 学 函数 的 关系 曲线 是 回 原点 的 扭 叶 型 。 分 析 
了 热 漏 系数 与 不 完全 回 热 因子 对 性 能 系数 、 致 冷 率 和 生态 当 
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Abstract: An irreversible quantum Stirling refrigeration model is established in which working substance ls a single 
particle contained in a one-dimensional infinite potential well, and the system consists of countless replicas. The 
occupation probabilities in each eigenstates of the particle are expressed by the thermal equilibrium Gibbs 
distributions. The cycle is made of two iso-energy branches and two isochoric branches and with heat leakage 
between heat reservoirs and imperfect regeneration. Analytical expressions of COP, cooling rate, entropy generation 
rate and ecological function are derived. The curve between cooling rate and ecological function is a loop-shaped 
one. Effects of heat leakage and imperfect regeneration on COP, cooling rate and ecological function are analyzed. 
Keywords: finite time thermodynamics; one-dimensional infinite potential well; quantum refrigeration cycle; 
ecological performance optimization 
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0 引言 


在 微 纳米 尺度 HL 了 、 低 温 或 高 密度 条 件 下 ， 需 要 P5。 一 些 学 者 也 研究 了 不 可 逆 因 素 对 量子 热力 循环 
考虑 热力 循环 工 质 的 量子 特性 。 应 用 有 限时 间 热 力 性 能 的 影响 。 金 晓 昌 等 2 考虑 高 、 低 温 热源 间 的 弱 
学 理论 和 方法 513， 可 以 建立 各 种 量子 循环 模型 06- 耦 联 作 用 ， 引 入 了 筋 通 热 漏 。Feldmann 和 
24], 研究 其 热力 性 质 与 性 能 参数 ， 优 化 其 运行 Koslofft271 引 入 了 内 摩擦 系数 来 描述 量子 非 绝 热 现 

间 。 实 际 热机 中 ， 存 在 着 热 阻 、 热 漏 、 摩 擦 、 涡 象 。 吴 锋 等 Pa 考虑 热 漏 、 热 阻 以 及 有 限 速 率 传 热 对 
流 、 惯 性 效应 及 非 平衡 的 影响 ， 循 环 为 不 可 逆 循 环 不 可 道 量子 斯 特 林 制冷 机 性 能 的 影响 ， 金 晓 昌 [9 从 火 


xl 


陶 济 学 角度 研究 了 量子 斯 特 林 制 冷 机 最 优 性 能 。 
在 文献 [3,10,14,16,23,28,30] 的 基础 上 ， 本 文 将 
建立 一 个 包含 高 低温 热源 间 的 热 漏 与 不 完全 回 热 的 
不 可 逆 斯 特 林 制 冷 机 循环 ， 其 工 质 为 无 数 个 囚禁 于 
一 维 无 限 深 势 阱 中 的 粒子 。 导 出 循环 的 性 能 系数 、 
致 冷 率 、 糯 产 率 、 生 态 学 函数 ， 分 析 这 些 重 要 的 性 
能 参数 与 热 漏 率 和 不 完全 回 热 因 子 之 间 的 关系 。 


1 系统 的 量子 力学 描述 


本 文 研究 的 对 象 是 一 个 因 禁 于 宽度 为 L 的 一 维 
无 限 深 势 阱 中 的 粒子 。 其 定 态 波 函 数 满足 : 


2 2 ei (D 
d“¥/dx +2mEYT/h =0 


根据 波 函 数 到 必须 满足 的 边界 条 件 WO)=0 和 
VD=0， 可 求 出 波 函 数 为 : 


YO)= Pap O) 
式 中 mw) 是 归 一 化 的 本 征 函数 ; 
四 (CO = V277sin(azxy 万 
系数 a 满足 归 一 化 条 件 : 
Pll = p,=1 (4) 


式 中 mr=o2， 为 概率 密度 ， 系 统 的 能 量 本 征 值 为 : 


E, =nxh? /((2mL) B) 
式 中 m 入 分 别 为 粒子 的 质量 与 量子 数 ， 哈 密 顿 量 
的 期 望 值 为 : 


E=WalB,=YpE, (0 
n=. n=] 
经 典 循环 中 ， 活 塞 在 气缸 中 做 往复 运动 假设 
势 阶 壁 像 经 典 循环 中 的 活塞 一 样 运动 ， 当 势 阶 宽 度 
变化 时 ， 系 统 的 波 函 数 与 能 级 都 会 随 着 工 变化 。 
定义 广义 力 为 


(7) 
Y=-aW /dady, 


式 中 ，y 是 与 广义 力 ,对 应 的 广义 坐标 。 因 此 ， 作 
用 在 势 阱 壁 的 力 可 以 表示 成 : 
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F=-dE/dL = plnah? / (mp)] (8) 
n=1 

2. 量子 斯 特 林 致 冷 循环 

研究 因 禁 于 一 维 无 限 深 势 阱 中 的 粒子 ， 为 简单 
起 见 ， 只 考虑 两 个 能 级 。 循 环 工 质 为 无 数 个 这 样 的 
粒子 ， 每 一 个 都 被 办 禁 在 一 个 一 维 无 限 深 势 阱 中 ， 
粒子 在 能 级 上 的 占有 几率 由 吉 布 斯 分 布 决 定 。 循 环 
由 两 个 等 能 量 过 程 与 两 个 等 势 阱 壁 的 回 热 过 程 组 
成 ， 分 别 与 经 典 循环 中 的 等 温 过 程 与 等 容 过 程 对 
应 。 考 虑 不 完全 回 热 与 高 低温 热源 之 间 的 热 漏 ， 系 
统 从 温度 为 7 的 低温 热源 吸收 热量 ， 放 出 热量 给 温 
度 为 Ty 的 高 温 热 源 ， 如 图 1 所 示 。 该 循环 是 一 个 不 
可 道 量子 斯 特 林 致 冷 循环 。 过 程 1-2 为 等 能 量 过 
程 ， 在 此 过 程 中 系统 的 总 能 量 保持 不 变 ， 由 此 可 


em 


E; = DIE t+ pb = 

pnEn t+ prEr = pbk + pybk,, (9) 
式 中 Pi,B (二 1,2,L;s=1,2) 分 别 表示 概率 密度 与 对 
应 能 级 的 能 量 本 征 值 ,i=L 表示 任意 态 ，n 为 量 
数 。 把 式 (4) 和 (5) 代 入 式 (9) 可 得 : 


(4-3p1)/ DF =(4-3pa)/T =(4-3py)/ Db 
(10) 


图 1 量子 斯 特 林 致 冷 循环 示意 图 . 

Figure 1. Sketch of a quantum Stirling refrigeration 
cycle. 

由 式 (8) 和 (10) 可 以 求 得 在 等 能 量 过 程 1-2 中 ， 作 用 

在 势 阱 壁 上 的 力 为 : 


11 
F,=Xh (4-3p)/ (mLD) EY 


等 能 量 过 程 中 系统 的 总 能 量 保 持 不 变 ， 系 统 从 低温 
热源 吸收 的 热量 全 部 用 来 对 外 做 功 ， 因 此 : 


% zh 1 aL 
=W,=|———(4-3p,)— 
CO， 12 | E Pu L 2) 


=X (4-3p)mL) In(L, /LL) 
同 理 ， 在 过 程 3-4 中 ， 系 统 释 放 给 高 温 热 源 的 
热量 与 外 界 对 系统 所 做 的 功 等 于 : 


入 
2 =|W,|= | 一 — 
L 


=A (4-3py)mD) In(L, /5) 
过 程 2-3 与 4-1 为 回 热 过 程 ， 理 想 回 热 过 程 
中 ， 在 回 热 器 中 交换 的 热量 分 别 为 : 


上 


1 1 
BP 


Axoh 1 天 
0 EP) EIp) (14) 
2? 1 万 
人 lM-3p)- -3p0)] ts) 
五 
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从 式 (14) 和 (15) 可 以 看 出 ，Q23=Q41, 这 说 明 ， 理 
想 量 子 斯 特 林 制 冷 机 可 以 实现 完全 回 热 。 假 设 高 低 
温 热 源 之 间 存 在 热 漏 ， 热 漏 率 由 下 式 决 定 : 


O = za (mL?) 0 


式 中 o 为 常数 。 则 每 一 个 循环 的 漏 热量 为 : 


2 = 0 7 (17) 
式 中 ?为 循环 周期 。 在 4-1 这 个 回 热 过 程 中 ， 如 果 是 
理想 回执 过程， 系统 会 在 状态 1 结束 回 热 过 程 ， 由 
于 内 部 不 可 道 性 ， 实 际 的 斯 特 林 制 冷 机 存在 回 热 损 
失 ， 因 为 非 理想 回 热 ， 系 统 在 状态 1” 结 束 回 热 过 
程 。 引 入 一 个 非 理想 回 热 因 子 j(y>0)， 则 非 理想 回 
热 过 程 中 的 损失 可 表示 为 : 

2 = HO (18) 

因此 ， 从 低温 热源 吸收 的 净 热 量 与 释放 给 高 温 
热源 的 净 热 量 分 别 可 以 表示 为 : 

2 = Qu +O,+0, (19) 


O; 二 QO 本 2 下 2 (20) 
粒子 在 茶 一 能 级 上 的 概率 密度 由 吉 布 斯 分 布 
为 : 


Pl=1/{1+exp[z hh /KT 2mL)]} CD 
P=1/{+exp[z hh /(x KT, 2mL )]} (02) 
Pi=1/{1+exp[z hh /(x kT 2mL )]} (23) 
P=1/{1+exp[z hi /KT,2mL)]} (24) 


3.， 循环 周期 
假设 势 阱 壁 以 速度 v(D) 运 动 ， 平 均 速度 为 久 可 


得 : 
万 一 
L,-h=| vvDd=vwa 03) 
等 能 量 过 程 1-2 与 3-4 所 需 的 时 间 为 : 
T=T4= (LL)/v (26) 
过 程 2-3 与 4-1 为 等 势 阱 壁 过 程 ， 在 此 过 程 中 


系统 内 能 变化 。 假 设 内 能 随时 间 的 变化 率 为 常数 ， 
因为 在 势 阱 宽度 不 变 时 ， 系 统 内 能 是 温度 的 单 值 函 
数 ， 故 有 : 


dT/dt = +M.; 
式 中 M 是 一 个 与 系统 能 量 无 关 的 常数 ， 它 只 取决 于 
回 热 材料 的 回 热 特性 。 正 负 号 分 别 表示 过 程 2- 
3(i=1) 与 过 程 4-1(i=2)， 为 简单 起 见 ， 假 设 WMA)=M2。 


因此 ， 回 热 过 程 2-3 与 4-1 所 需 的 时 间 可 分 别 表示 
为 : 
(28) 


t=(T -7.)/M 
t=(T, —T.)/M (29) 
整个 循环 的 周期 可 以 表示 为 : 


7 三 Th 十 723 十 934 十 2 三 


2(D -DL)/v+2(T, -7T,)/M 
4. 不 可 道 斯 特 林 制 冷 循环 一 些 重 
能 参数 及 分 析 


由 式 (12〉 和 (13) 可 得 循环 的 输入 功 为 : 


(30) 


要 性 


XA"h? a L 人 联 立 式 (19) 与 (30) 可 得 的 致 冷 率 为 ， 
mL Xx R={(4-3pi)Inx—a[2L (x—D/v+ 
式 中 x=LyL1， 表 示 势 阱 宽度 比 。 由 式 (19) 和 (31) 可 2(T5 —T)/M]-A(4-3p1) 
得 该 不 可 北 量 子 斯 特 林 制冷 机 的 性 能 系数 (COP) (4—3p3)/x’]}/{2L (x-1)/y 
为 : +2(7 -7T.)/M} 

cop={(4—3pi)lInx—ar—ul(4—3p) 
—(4—3p3)/x ]}/{Inx{(4—3p,,) (32) 
—(4—-3pa)/x ]} 


_1z 
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W = 


(33) 


炳 产 率 为 


-4- 3pu) lnx— CT 一 AL[(4 一 3p11)— (4— 3pa) /x De 3p31) 


kmL? 
eA pe a i 


Ax’h? 


i ~——(1— pai)] 


kmL kmL 
2T,x i D2)]+3[ 一 一 人 2T,x 二 Di) 一 3[ 一 一 六 元 2 万? 2T,x 


xh? 


类 输出 率 为 : 


1 ,TT ,x 

zt{(7 1) mL {(4-3pi) Inx-or— Hl(4-3p) 
L 

1 


Mk 古 前 yo 
二 
(35) 
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(4-3ps)/ + (了 Do Ct 3pa)Inx+ar+H(4-3pn)—(4-3p3a)/x 1]}} 


生态 学 目标 函数 B133] 为 : 


4 1 TT 元 大 T 
E=—-oNh=—{(2 -Ds{(4-3pu)nx -oar -AH(4-3p1)-(4-3pu)/x t(D) 
T 7 7 mL T, 


xh 1 二 zh? 


] 
7 {-—(4-3pa)lnx+art+nul(4—-3p1)—(4-3p3)/x])} 7 {—{(4—3pi)Inx— 
mL (I 


1, 1 
or — Hl(4—3p1)—(4 3p)/x +i {4-3p nxtortA(4-3pn)-(4-3p) (36) 
H 


kmL? 1 kmL? 
a J ne 


1 km 
十 a7 — pa)]+ i 十 


1 kmL 1 
2 ( pi) 3 ln Pa + 


SE = EmL /1(z272) 为 无 量 纲 生 态 学 目标 函 率 、COP 和 生态 学 目标 函数 都 是 势 阱 宽度 比 x 的 函 
数 。 图 2 给 出 了 非 理 想 回 热 因子 w 取 不 同 的 值 时 ， 
性 能 系数 随 x 的 变化 关系 。 计 算 中 取 


数 。 由 式 (27), (28) 和 (31) 可 知 ， 对 于 给 定 的 Ts, 


To LM VQ WC = 如 /(m 忆 ), 制 冷 机 的 制冷 
re A SH OO M=15, M =15,2/(L,v)=1, 


chinaxiv:201801.00075V1 


0 =0.01,C/k=300。 由 图 可 知 ， 制 冷 机 的 性 能 
系数 随 势 阱 宽度 比 x 的 增加 而 减 小 。 随 着 非 理想 回 
热 因 子 w 的 增加 ，COP 减 小 。 图 3 给 出 了 /=0.01， 

a 取 不 同 的 值 时 ，COP 随 x 的 变化 关系 。 计 算 中 参 
数 取 值 与 图 2 相同 。 从 图 可 知 ，COP 随 着 a 的 增加 


EN ,Re 
Oe Ul = 


24 
3 
Xx 
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而 减 小 。 图 4 给 出 了 制冷 率 随 x 的 变化 关系 ， 计 算 
中 c=0.1， ”kX 取 不 同 的 数值 ， 其 他 参数 的 取 值 与 图 
2 相同 。 从 图 可 知 ， 制 冷 率 随 x 的 变化 关系 为 类 抛 
物 线 关 系 ， 由 极 值 条 件 OR /Ox = 0 得 : 


,mL 1 
[(4—3p3)+3p3 mel 


(37) 


-4 3p Inx or-— Hl(4-3p1) 三 (4 3p31)]}=0 


式 (36) 给 出 了 制冷 率 与 势 阱 宽度 比 x 间 优 化 关 
系 。 图 5 给 出 了 /=0.01，a 取 不 同 的 数值 时 ， 制 冷 
率 随 x 的 变化 关系 ， 制 冷 率 随 着 热 漏 的 增 大 而 降 
低 。 


图 6 给 出 了 无 量 纲 生态 学 目标 函数 与 制冷 率 
之 间 的 优化 关系 (BZ*-R)。 从 图 可 知 ，E”-R 优化 关系 
曲线 为 扭 叶 型 ， 存 在 着 一 个 生态 学 目标 函数 的 极 大 
值 点 与 一 个 制冷 率 的 极 大 值 点 ， 且 每 一 个 生态 学 目 
标 函 数值 ( 除 极 大 值 点 外 ) 对 应 两 个 制冷 率 取 值 ， 显 
然 要 取 对 应 的 制冷 率 较 大 的 状态 点 。 热 漏 不 改变 曲 
线 的 类 型 ， 但 随 着 热 漏 的 增加 ， 生 态 学 目标 函数 值 
和 制冷 率 都 变 小 ， 在 生态 学 目标 函数 取 值 较 小 时 ， 
热 漏 对 其 影响 较 小 ， 在 生态 学 目标 函数 取 值 较 大 
时 ， 热 漏 对 其 影响 较 大 。 图 7 给 出 了 对 E*-R 关系 
曲线 的 影响 ， 比 较 图 6 和 7 可 知 ， 热 漏 对 E*-R 关系 
的 影响 比 非 理想 回 热 因子 的 影响 大 ， 即 热 漏 是 比 非 
理想 回 热 更 重要 的 不 可 逆 因 素 。 
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图 2 /对 COP 与 势 阱 宽度 比 的 关系 的 影响 . 
Figure 2. Effects of 4/ on COP with the ratio of potential 


well width. 
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图 3 a 对 COP 与 势 阱 宽度 比 的 关系 的 影响 . 
Figure 3. Effects of ac on COP with the ratio of 


potential well width. 
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4. /对 致 冷 率 与 势 阱 宽度 比 关系 的 影响 . 


Figure 4. Effects of nu on cooling rate(R) with the ratio 


ofpotential well width 
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本 文 建立 了 以 无 数 个 无 限 深 势 阱 中 的 粒子 为 工 30 
质 的 不 可 逆 量 子 斯 特 林 致 冷 循 环 模型 。 基 于 薛 定 雇 | 1 si 


| 冷 机 的 性 3.a=0.05h=01 


20r 


方程 的 求解 分 析 了 该 不 可 逆 量 子 斯 特 林 
能 ， 给 出 了 性 能 系数 、 致 冷 率 、 炉 产 率 和 生态 学 目 
标 函 数 ， 分 析 了 热 漏 系 数 与 不 完全 回 热 因子 对 性 能 


it 


系数 、 致 冷 率 以 及 生态 学 函数 的 影响 。 制 冷 率 与 生 | ] 
1.4 5 上 ]] 
12r 1.0=0,u=0.01 二 
1 2.0=0.1,u=0.01 0 
1 3.0.=0.3,4=0.01 | 0 01 02 03 04 05 06 07 08 0.9 1 
4.0=0.5,u=0.01 民 


图 7 kh 对 致 冷 率 与 无 量 纲 生态 学 函数 关系 的 影 
啊 . 


Figure 7. Effects of 1 on cooling rate versus 


dimensionless ecological objective function. 
《制冷 率 与 生态 学 函数 的 关系 曲线 是 回 原点 的 扭 叶 
图 5 oa 对 致 冷 率 与 势 阱 宽度 比 的 关系 的 影响 ， 型 。 存 在 者 一 个 生态 学 目标 函数 极 大 值 与 一 个 制冷 
率 的 极 大 值 。 这 两 个 极 大 值 之 间 区 域 是 该 制冷 机 的 
ratio of potential well width . 最 优 运行 区 间 。 本 文 结 果 有 助 于 加 深 对 以 无 限 深 势 
ao 阱 中 粒子 为 工 质 的 不 可 逆 量 子 斯 特 林 制冷 循环 性 能 
ee 的 理解 。 
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Figure 9. Effects of o on cooling rate versus the 
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